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Overview & Motivation 



• Overview 

- Presentation of computational and experimental 
results from a recently developed rectangular wing 
aeroservoelastic modeling tool 

• Motivation 

- Compare tool to independently published work^ 

- To support rapid investigation of aeroservoelastic 
phenomena in a medium-fidelity tool 

• Also novel sensors such as fiber optics 

- Provide a rapid aeroservoelastic design platform 
which can serve students of aeroservoe asticity 

^SConyers, H. J., Dowell, E. H., and Hall, K. C., Aeroservoelastic Studies of a Rectangular Wing with a Hole: 
Correlation of Theory and Experiment," 2010 Aerospace Systems Conference 
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Background 


• In previous work^ tool used to model a clamped 
wing structure with two control surfaces and fiber 
optic sensor feedback used for flutter suppression 


Computational Wing Model 


Active Flutter Suppression 


Wing Root 


Three fibers with Bragg gratings spaced 
every 14 inch 



X, m 
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Aeroservoelastic Tool Overview 


Tool for Rectangular Wing Aeroservoelastic Design in MATLAB 


Tool allows the user 
to quickly move from 
inputs like aspect 
ratio, control surface 
count, and half span 
to a linear time 
invariant state space 
model which can be 
used for control 

- A few seconds of real 
time computation 

- Most important 
structural and 
aerodynamic properties 
are parametric 


Geometry Pararneters 

Aspect Ratio 
Half Span 

Spar/Rib/Skin dimensions 

Control Surtece 
NumbedSize/Layout 

All Materials 

Oiscretizaton of Model (x,y ) 
Aerodynamic Parameters 


Reduced Frequency Range 

Lag constants 

N u m be r of a e rody nami c 
lag states 


State Space Parameters 


Actuator T ran star F u notion 



12/17/2014 


AIAA Modeling and Simulation 
Technologies Conference, Jan 5-9, 2015 


5 



Graphical Path of Verification and 

Validation of Tool 


Compare to State Space 
Models 



Generalized 
Aero Forces 


Compare to Wind Tunnel 
Test 
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Direct 

Comparison 


Rational 
Function 
Approximation 







Beam Model Verification 


• Beams used to model wing structure 

- FEM with 30 elements compared to theory 
show good matches in bending and torsion 



Cantilever Theory versus FEM Beam Model: Deflection and Twist 



Graphical Path of Verification and 

Validation of Tool 








^SConyers, H. J., Dowell, E. H., and Hall, K. C., Aeroservoelastic Studies of a Rectangular Wing with a Hole: 
Correlation of Theory and Experiment," 2010 Aerospace Systems Conference 


Plate FEM Validation 



• Ground Vibration Test (GVT) on a 
article used for validation of plate 
FEM 

• Plate FEM Discretized with 16x16 
12 DOF isotropic plate elements 

• Experiment shows good 
correlation with ANSYS and tool 


GVT on Article (with a hole 
for a different test) 



Computational Model 

I = 304.8mm 







16x16 elements 

3^ 


■o ^ 

th = 1.588m 

fD y 



Q- ! 




Experimental Data^ 



ANSYS 
Frequencies, Hz 

Tool FEM 
Frequencies, Hz 

Conyers et al. GVT, Hz 

Mode # 1 

3.99 

3.99 

4.13 

Mode # 2 

16.96 

16.97 

17.24 

Mode # 3 

24.86 

24.89 

24.38 

Mode # 4 

55.33 

55.40 

54.25 

Mode # 5 

69.84 

69.92 

69.00 
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152.4mm 



Graphical Path of Verification and 

Validation of Tool 


Compare to State Space 
Models 



Compare to Wind Tunnel 
Test 
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RFA Verification 


Generalized 
aerodynamic forces 
(GAF) computed for 
plate 

Roger's rational 
function approximation 
(RFA) used to fit GAF 
coefficients a> 

- 4 lag states 

Least squares error for i 
bending and twist | 
coefficients 


Generalized Aerodynamic Force 

Q(ik) = Zf'^ Ditk)-^ApWc,p, 

Rational Function Approximation of GAF 


Q(s) — Aq -\- sAi + s^A2 + 


I 

1=1 


+ ft 


^2 + 1 


Comparison of GAF and RFA Curve Fits 
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^ First bending GAF 
First bending RFA 

k = 0.8 


k = 0.6 


k = 0.4 


k = 0.2 


b) 


16 r 


14 ■ 


12 ■ 


10 - 


8 

6 ■ 
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0 - 
-35 



-34 -33 -32 -31 -30 -29 -28 


A- = 0.4 


A = 0.2 * 


First torsion GAF 
First torsion RFA 


Real part 


Real part 
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V-g Analysis using RFA 



• The test plate article flutter speed was predicted to be 19.9 
m/s 

— traditional bending/torsion flutter mode 


V-g Analysis on Computational Plate Article 
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V-f Analysis using RFA 

• The test plate article flutter frequency was predicted to 
be 10.9 rad/s 

— Torsional mode shifts closer to bending mode 
— Characteristic of a one side clamped plate flutter mode 



V-f Analysis on Computational Plate Article 
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Graphical Path of Verification and 

Validation of Tool 


Compare to State Space 
Models 



Generalized 
Aero Forces 
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^SConyers, H. J., Dowell, E. H., and Hall, K. C., Aeroservoelastic Studies of a Rectangular Wing with a Hole: 
Correlation of Theory and Experiment," 2010 Aerospace Systems Conference 


Flutter Validation Experimental Stud 





A wind tunnel investigation was 
completed at Duke University in 
previous work 

— Tool flutter speed shows good 
correlation with Conyers et al.'s flutter 
code 

• Differences may be due to use of more 
aero panels in the tool 

— Wind tunnel results were comparably 
close 



Configuration for Wind Tunnel Test at Duke 
(different article shown than flat plate) 


Experimental Data^ 


Conyers et al. 
Flutter Code^ 

Tool Flutter Code 

Conyers et al. Wind 
Tunnel Results^ 

Flutter speed , m/s 

20.8 

19.9 

20.05 

Flutter frequency, Hz 

10.3 

10.9 

11.50 
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Graphical Path of Verification and 

Validation of Tool 



Compare to Wind Tunnel 
Test 
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State Space Model Verification 

• We verify that the state space models 
correlate with what was predicted from the V- 
g and V-f analyses 
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State Space Model Architecture 


Components of state space models 

- FEM mass, stiffness, damping and modal matrices 

- Rational function approximation coefficients 

- Actuator dynamic models 

- Flight condition 

State Space Model Architecture 
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Analytical Model with Control Surface 

• Verification of state space models is completed for a wing 
model with 

- internal aluminum beam spar and rib structure 

- aluminum skin 

- a control surface and a leading edge accelerometer 
Analytical Model with One Control Surface and a Leading Edge Accelerometer 



3.0 0-2 X, m 
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Actuator Dynamics 

• Actuators are modeled as order transfer functions 

- 1^* order command lag 

- 2"'' order actuator dynamics 

Actuator Model with and without command lag 
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Analytical Wing Mode Shapes 

• Mass normalized mode shapes are computed with high 
torsional spring stiffness in connected control surfaces 

• Control modes are computed with low torsional spring 
stiffness and a prescribed 1 deg. rotation boundary condition 


Analytical Wing Modal Analysis 

a) First bending, frequency = 16 rad/s b) First torsion, frequency = 43.4 rad/s 



c) Second bending, frequency = 85.6 rad/s d) Second torsion, frequency = 184.6 rad/s 




State Space Model Verification 

• We verify that the state space models 
correlate with what was predicted from the V- 
g and V-f analyses 
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V-g Analysis with RFA 


• V-g analysis of wing shows a traditional 
bending/torsion flutter mode appearing at 76.5 
m/s 

V-g Analysis of Analytical Wing Model 
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Wing Model Pole Migration 


• The bending mode becomes more stable 

• The torsion mode becomes neutrally stable at 76.5 m/s 

• Flutter speed is the same as predicted in the V-g analysis 


Pole Migration of State Space Model from 20 - 78 m/s 
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State Space Model Verification 

• We verify that the state space models 
correlate with what was predicted from the V- 
g and V-f analyses 
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V-f Analysis 



• Frequency analysis shows the flutter 
frequency at 28 rad/s 


V-f Analysis of Analytical Wing Model 



Bode Plot of State Space Model 

• At speed below flutter speed, amplitudes of two distinct modes visible 

• At flutter speed only flutter mode is visible 

• Frequency is the same as predicted from the V-f analysis 


Bode Plot of Surface to Leading Edge Accelerometer 




State Space Model Verification 



• We verify that the state space models 
correlate with what was predicted from the V- 
g and V-f analyses 
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Impulse to State Space Model 


• Flutter is apparent in model designed past flutter speed 
- Divergent oscillatory 

• Model at lower speed is damped after impulse 


25 
20 1 - 
15 


V = 77m/s 

V = 40 m/s 



Effects of higher 
modes are visible 
before flutter 
dominates 


~4.5 peaks per 
second 
~28 rad/s 
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State Space Model Verification 



• We verify that the state space models 
correlate with what was predicted from the V- 
g and V-f analyses 
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Velocity, m/s 


1-COS Gust Model 

Gust inputs to structure are designed with gust 
modes and 1-cos gust input structure 



1-cos Gust Input Model 



Gust mode approximation 

2k 

( -V 

c.p. 


3 wash ~ I ^ - (^c.p. ^gust 


) 


Gust velocity 
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w, 


g.max 


(1 — cos(Gt)) 


Gust acceleration 


Wa(0 = 


W, 


g,max 


Time, s 


2 

Gust frequency 


sin(Gt) G 
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Gust Input to State Space Model 

• The response of wing to 1-cos gust is expected 

— Low frequency gust response and high frequency 
oscillations from flutter are seen to be superimposed 
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Conclusions 



• Several first step verification and validation 
studies were presented for a new 
aeroservoelastic tool 

• More verification and validation is needed to 
assess the state space models including 

— An experimental flutter test and active flutter 
suppression 

• This work further supports independent flutter 
analysis conducted by Dr. Conyers in his 
dissertation 
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Future Work 



• Improvements will be made to include rigid 
body modes in the tool 

• Input structure will be made more user 
friendly 

• Would like to look into transitioning to use as 
an open tool for students 
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Questions? 
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